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Abstract 

Wc have simulated a granular medium using the Discrete Element Method to study the influence of particle 
shape on the behaviour of the whole assembly of 5,000 grains. By comparing two shape groups: polygons 
and clumps of discs, we observed the effects of grain shape variations on both macro and microscopic levels 
and investigate the causes. 
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Abstract 

We performed a scries of numerical vertical compression tests on assemblies of 2D granular material using 
a Discrete Element code and studied the results with regard to the grain shape. The samples consist of 
5,000 grains made from either 3 overlapping discs (clumps - grains with concavities) or six-edged polygons 
(convex grains). These two grain type have similar external envelopes, which is a function of a geometrical 
parameter a. 

In this paper, the numerical procedure applied is briefly presented followed by the description of the 
granular model used. Observations and mechanical analysis of dense and loose granular assemblies under 
isotropic loading arc made. The mechanical response of our numerical granular samples is studied in the 
framework of the classical vertical compression test with constant lateral stress (biaxial test). The compar- 
ison of macroscopic responses of dense and loose samples with various grain shapes shows that when a is 
considered a concavity parameter, it is therefore a relevant variable for increasing mechanical performances 
of dense samples. When a is considered an envelope deviation from perfect sphericity, it can control me- 
chanical performances for large strains. Finally, we present some remarks concerning the kinematics of the 
deformed samples: while some polygon samples subjected to a vertical compression present large damage 
zones (any polygon shape), dense samples made of clumps always exhibit thin reflecting shear bands. 
This paper was written as part of a CEGEO research project^] 
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1. Introduction 

A typical numerical approach to discrete element modeling of granular materials is to use simple shapes 
of particles (discs in 2D [l[ or spheres in 3D Q). Although the computation time is short using this method, 
these models cannot reflect some of the more complex aspects of real granular media behaviour, such as 
high shear resistance or high volumetric changes |3j. In order to model these mechanisms properly, physical 
phenomena (resistance to inter-granular rolling [J, [H, Q) or other grain shapes (sphere aggregates p} or 
polyhedral grains @) must be used. The influence of grain shape is not yet fully understood. In this article, 
we will present our findings concerning the influence of grain shape (grain concavity in particular) on the 
mechanical behaviour of granular assemblies. We compared two groups of grains - convex irregular polygons 
and non-convex aggregates of three overlapping discs. 
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2. Granular Model 



The granular model used consists of 5,000 polydisperse 2D frictional particles. Two kinds of grain shape 
are used: convex irregular polygons with six edges and non-convex particles made of aggregates of three 
overlapping discs called clumps. These two shapes were chosen because of the similarity of their global 
contour (polygonal grains can be seen as a polygonal envelope of clumps made of three discs). As shown 
in Fig. ??, particle shape is defined by a parameter a = ^jp, where R\ denotes the ex-circle radius of the 
particle and AR is the difference between the ex-circle and the in-circle radii, 0, [13]. The in-circle must be 
fully contained in the particle. For non-convex clumps, a ranges from (circle) to 0.5. For convex polygonal 
grains, a ranges from 1 — ^ ~ 0.13 (regular hexagons) to 0.5 (equilateral triangles). Some of the shapes 
used are presented at the bottom of Fig. ??. For each chosen a, granular samples are made of polydisperse 
particles: the polydispersity of grains is determined by the radii of the grain ex-circle. In each sample, the 
chosen radii R\ are such that the areas of the ex-circles are equally distributed between S m = ir(R m ) 2 and 
S M = k(Rm) 2 = ^(3i? r „) 2 . 



3. Discrete Element Method 



Two-dimensional numerical simulations were carried out using the Discrete Element Method according to 
the principles of Molecular Dynamics (MD) Two codes were used: PFC 2D by ITASCA [H] for clump 
simulation and a code capable of dealing with polygonal particles that was developed at the laboratory. 
Both codes use the same contact laws for contact forces computations [l3|: grains interact in their contact 
points with a linear elastic law and Coulomb friction. The normal contact force /„ is related to the normal 
interpenetration (or overlap) h of the contact 

} n = k n -h , (1) 

as f n vanishes if contact disappears, i.e. h = 0. The tangential component f t of the contact force is 
proportional to the tangential elastic relative displacement, according to a stiffness coefficient k t . The 
Coulomb condition \f t \ < fi f n requires an incremental evaluation of f t in every time step, which leads to 
some amount of slip each time one of the equalities ft = ifJ-fn is imposed (fi corresponds to the contact 
friction coefficient). A normal viscous component as opposed to the relative normal motion of any pair 
of grains in contact is also added to the elastic force /„. Such a term is often introduced to facilitate the 
mechanical equilibrium approach [l4| . In case of frictional assemblies under quasi-static loading, the influence 
of this viscous force (which is proportional to the normal relative velocity, using a damping coefficient g n ) is 
not significant [l5| (elastic energy is mainly dissipated by Coulomb friction) . Finally, the motion of grains 
is calculated by solving Newton's equations using either a leap-frog (in PFC 2£> ) or third-order predictor- 
corrector discretisation scheme [l6| (in the in- house software). This constitutes the only known difference 
between the two codes. 

The principles of disc contact detection are well known [l^, [13] , and contact detection for clumps was 
solved in the same way: contact occurs at a point, the normal force value /„ is calculated with cq. (H} 
and its direction connects the centers of discs in contact, Fig. ??. Contact detection and contact force 
calculations between polygons do not use classical methods based on the area overlap between polygons [3, 
SHI Thc shadow overlap technique proposed by J.J. Moreau [22| . which was originally applied within 



the Contact Dynamic approach [231 ] for convex polygonal particles, was used. In our study this technique 
was adapted to the MD approach. Three types of geometrical contact can exist between polygons: Corner- 
to-Corner, Corner-to-Edge (CE) and Edge-to-Edge (EE). Corner-to-Corner contacts are geometrically (or 
mathematically) realistic but never occur in our simulations because of the numerical rounding errors. When 
dealing with (EE) contact, contact detection involves two contact points and their associated overlaps h, 
Fig. ??. This is the main difference compared to thc classical method (area overlap calculations) where only 
one contact is considered between the edges. 

Finally, we may be interested in the main contact law parameters: the normal and tangential stiffness, k n 
and k t , and the friction coefficient fi. Assuming that samples would first be loaded with a 2D isotropic stress 



the normal stiffness of contact k n was calculated according to the dimensionlcss 2D stiffness 
kn/co 15, 24, 25, . k expresses the mean level of contact deformation, l/n = h/ (2i?), 



a = 10 kN/m, 
parameter k = 

where (R) is the mean particle radius. In our simulations, k was arbitrarily set to 1,000. As a comparison, 
a sample made of glass beams under isotropic loading of 100 kPa reaches k = 3, 000. The tangential stiffness 
kt can be expressed as a fraction of the normal stiffness, k = k t /k ni k > 0. k > 1 may exhibit specific 
behaviour where Poisson coefficient of grain assemblies become negative 27, 28, 2jJ 3^]. Running several 
numerical simulations with various k, < k < 1, j25| have shown that if 0.5 < k < 1, the macroscopic 
behaviour remains similar. Thus we arbitrarily set k to 1. 



4. Sample Preparation - Isotropic Compression 

Granular samples of 5, 000 grains are prepared in three steps: preparations start with a random spatial 
distribution of particle position inside a square made of four rigid walls. Secondly, the particles expand 
slowly until oo = 0.5 kN/m is reached. Finally, samples are isotropically loaded by wall displacement up to 
o"o = 10 kN/m. To obtain samples with different compacities, we may use various values of the inter-granular 
friction coefficient fi during the preparation [3lj ]. When /i is set to zero, samples isotropically loaded up to 
o"o = 10 kN/m are dense and the compacity is maximal. When a strictly positive value of /i is used instead, 
samples become looser and compacity decreases. In our study dense samples were prepared with [i = and 
the loose ones with \i equal 0.5. 16 different samples (4 dense, 4 loose made of clumps and 4 dense, 4 loose 
made of polygons) for each a value were prepared_|. Dense samples will be written as Crj or Prj respectively 
for Clumps and Polygons. Loose samples will be denoted Cl or Pl . A subscript can be added. It then 
corresponds to the decimal part of the shape number a. As an example, fragments of two dense samples 
with a = 0.30, Cd-30 and Pd-30 , are displayed in figure ??. 

For both clumps and polygons, contact between particles can occur at more than one contact point. 
There are four contact possibilities for clumps: single contact Fig. ??(i), double contact involving three 
discs Fig. ??(ii), double contact involving four discs Fig.??(iii) and triple contact Fig. ??(iv). By analogy 
with polygon contacts (Edge-to-Edge or Corner-to-Edge, Fig.??), all these contacts between clumps can be 
merged into two groups, (CE) and (EE). In the (CE) group, grains involved in a single contact (i) may 
rotate without sliding. Double contact (ii) allows rotation with sliding and eventually friction. Rotation 
and sliding of polygons meeting at a single contact point are not correlated. Group (EE) contacts block the 
rotation of the grains. In the case of rotation of grains, contacts of this type would be lost. Therefore, the 
shape of the grains can be regarded as macro roughness. 

For samples subjected to isotropic loading, we focus on two internal parameters that mainly determine 
the mechanical behaviour: the compacity £ and the coordination number z* . Figure ?? shows the evolution 
of £ with a for dense and for loose samples. For Cpj samples, £ evolution is bell-shaped and maximum 
compacity is reached for 0.2 < a < 0.3. Similar observations were made by This also seems to be the 

case for Cl samples although the amplitude of the bell-shaped curve appears to be lower. The geometrical 
origin of these results deals with the grain shape (concavity and grains envelope) and imbrication and the 
interlocking between grains, but appears to be complex to establish. While we might think that when 
two grains are in contact with a single point of contact (contact (CE)-type (i), Fig. ??), this would tend 
to increase the local porosity and thus reduce the overall compacity, this does not seem to be the case: 
figure ?? shows clearly that for Crj samples, the proportion of contact (CE)-type (i) is at its maximum for 
0.2 < a < 0.3 in the case of dense samples and decreases linearly with a in the case of Cl samples. This 
completely goes against the £ trend, even if it is clearly established that when £ is high (Cl> samples), the 
proportion of (CE)-type (i) contacts is lower than in the case where the compacity is low (Cl samples). 

For samples made of polygons, the dependence of £ on a is not clear either. For Ppj samples, £ is almost 
constant for 0.2 < a < 0.3. If we exclude the Pd-40 sample, which behaves in a peculiar way, we can notice 
that for samples made of grains with regular shapes, corresponding to a = 0.13 and a = 0.5, the compacity £ 



1 All the analyses presented in this article were carried out on mean results calculated over 4 samples of each density and 
each a. Associated Standard Deviations will always be given, even if they are too small to be significant 



is smaller than for samples made of grains which have shape irregularities (0.13 < a < 0.5). The relationship 
between £ and the percentage of (CE)-type (i) contact is again impossible to establish. Furthermore, it is 
worth observing that when 0.2 < a < 0.3, Crj and Prj samples show very similar compacity. Here, grain 
envelopes seems to be a cleverer interpretation of a parameter than grain imbrication, which is only relevant 
for clumps. Finally, we should note that the angle of friction used during the preparation does not seem to 
have a major influence on the compacity of samples made of polygons when 0.2 < a < 0.3. 

For samples subjected to isotropic loading, we studied the coordination number z* corresponding to 
the mean number of contacts per grain. Here only grains that have two or more compression forces, and 
therefore take part in the load transfer, were considered. For samples made of frictionless perfectly rigid 
discs, z* = 4 [32j. Because k is not infinite in our study, in the samples made of frictionless circular particles 
(Cpj.00 ), z* = 4.093 ± 0.005 is greater but still very close to the reference value, z* is evaluated for clumps 
and polygons and both dense and loose samples. The dependence of z* according to a is shown in Fig. ??. 
For Cl samples, z* increases linearly with a, like Pl samples, but for Cl , it can be directly correlated to 
the percentage of (CE)-type (i) contacts which decrease with a, and then increase z* . For dense samples, 
the percentage of (CE)-type (i) contacts did not vary too much, z* is constant for Crj and Ppj samples. 

5. Macromechanical Response of Granular Assembly Loaded in Vertical Compression Test 

The samples were tested in a 2D strain controlled vertical compression test, also called biaxial test. 
Vertical stress o\ was applied by increasing the compressive vertical strain e\ while lateral 03 remained 

constant. The loading velocities were chosen according to the dimensionless inertial number I = iiy^^- 



[261 ] where i\ denotes the strain rate and (m) the typical mass of a grain. It describes the level of dynamic 
effects in the sample. For quasi-static evolutions, the value of / should be low. / value was set to 5 • 10~ 5 
for clumps and for polygons samples, regardless of the code used. During the vertical compression in both 
dense and loose samples, the same value of contact friction coefficient [i = 0.5 was used. The mechanical 
responses of the samples are plotted on r\ vs. E\ charts and shown in figure ??. r\ ~ t/s, t = (<ti — o- 3 )/2 
is half of the deviator stress and s = (<ti + oj,)j2 is the mean stress. Extracted from r\-£\ curves shown 
in figure ??, friction angles at the peak <f> p and at the threshold (f>t are given figure ??. Average dilatancy 
angles extracted from Fig. ?? are presented in the figure ??. 

For dense samples made of clumps, Crj , we can observe in figure ?? that the macroscopic shear resistance 
increases with a. Although Crj-10 implies grains with a small a, the mechanical response of the sample 
exhibits remarkable increase of the maximum deviator in comparison to disc samples Cpj.OO where rotations 
of particles are not potentially disturbed by the grain shape. For Cl samples, there is no peak value of the 
friction angle, 4> p = </>t, figure ??. Thus, for all clump samples, both peak <p p and threshold 4>t friction angle 
values increase along with a. <pt increases proportionally with a while the increase of <f> p is nonlinear and 
seems to be asymptotic for a > 0.4. 

For Ppj samples, <j) p values slightly decrease linearly along a, while 4>t values increases, Fig. ??. 

For Pl samples, we can notice in figure ?? that the macroscopic curves show typical behaviour similar 
to that of loose samples when a < 0.3 and a typical behaviour characteristic of dense samples for a > 0.3. 
This kind of behaviour deals with the values of the initial compacity of the samples shown in figure ?? where 
we can observe that £ values for Ppj and Pl samples are very close when a < 0.3. If we focus only on <f> t 
for Pl samples, we can observe (Fig. ??) an increase of the friction angle with a, except for samples made 
of triangles, a = 0.5, which always behave in a specific wajH. In conclusion, it can be noted that adding 
some shape irregularity by increasing a always leads to an increase of the macroscopic angle of friction in 
the critical state. This is the case for clumps and polygons with a constant microscopic friction angle fi. 

This influence of grain geometry is in line with a previous study by Salot et al. 0. Lastly, as we can 
see in figure ??, a does not explicitly influence Young's modulus. E is linked to the rigidity matrix and 
therefore to z* [33J, which is constant for dense samples (Fig. ??). 




2 Note that triangle is the only shape with 3 edges 



Similarly, particle concavity does not particularly influence the average dilatancy angle tp values (sin i\) = 
del-del ) °^ dense and loose clump samples (Fig. ??). On the other hand, ip is lower for polygons with higher 
values of a (closer to triangular shape) than for those more similar to hexagons. 

In figure ??, volumetric changes in some samples are illustrated. For both dense clump and polygon 
samples, Fig. ?? and ??, the volumjl increases mainly during vertical compression, after a small contraction 
due to the stiffness of the contacts 15j. The volumetric increase for Crj samples is quite similar from one 
a to another (Fig. ??). On the other hand, a clearly influences the volumetric change of Pd samples but 
this influence seems erratic. Nevertheless, Pd samples show greater total dilatancy than Cd samples. 
Cl samples (Fig. ??) behave like loose sands and contract all throughout the compression test. It is more 
complex for Pl samples (Fig. ??) for which the setting-up process remains problematic for some values of 
a. 



6. Micromechanical Analysis 

From the macroscopic results exposed in the previous section, two main observations can be established: 
for dense samples made of clumps, the evolution of the angle of friction at the stress peak (f> p varies signifi- 
cantly with a. The geometrical imbrication between grains in contact, which depend on a, may be one of 
the micro-mechanical origins of these results. Secondly, for all samples, dense or loose, made of clumps or 
polygons, it was established that the angle of friction <fi t at the end of biaxial tests increases with a and 
is independent of the initial state. This result is proof of the role of the grain envelope in the mechanical 
behaviour rather than some inter-granular imbrication considerations. In this section we will try to gather 
evidence for these proposals. 

6.1. Contact proportion evolution for dense samples made of clumps 

Single contacts are mainly involved in all the grains samples tested, Fig. ?? {(CE) contacts for polygons 
or (CE)-type (i) contacts for clumps - Fig. ??). During loading, it can be noted that the percentage of single 
contacts increases. Because dense samples were prepared with no friction, compacity of each assemblies of 
grains subjected to isotropic loading is constant and always maximal. As a consequence, even if samples can 
exhibit slight contractancy (related to dimensionless contact stiffness n, [Hj]), the total number of contacts 
in each sample during a vertical compression systematically decreases. Therefore, in order to compare 
different contact type proportions in different phases of the test, we suggest balancing the decrease in the 
total number of contacts using the coefficient ui* = N* IN* . N* and N* represent the total number of 
neighbouring contact^ at the given vertical strains £\ (e a or £&) in the samples. Here, we suggest focusing 
on the evolution of two new contact groups for clumps: 

• (SC) single contact between grains, known as simple contact, 

• (CC) multiple contacts between grains, hereafter called complex contacts. 
These two new groups can be examined in figure ??. 

We observed the evolution of clump contact numbers of each group between two successive stages by 
normalising this evolution with uj* . We thus defined a new variable A = a;* • N 6a /N £b , where N ea and N eb 
denote the number of (SC) or (CC) on two different levels e a and Eb of the vertical strain e\. In figure ??, 
A evolution is calculated between e a = £j (isotropic state) and Eb = £ P (maximum stress deviator). We can 
observe that for a = 0.10, A is smaller than 1 for ( CC) and greater than 1 for (SC). This can be regarded as 
a transformation of (CC) into (SC) between these two stages. If all the complex contacts transformed into 
simple, graphic points were at an equal distance from 1,1 — ^(QC) ~ ^(SC) ~ ^ If 1 — ^(CC) ^ ^(SC) ~ 
it means that some complex contacts transform into simple contacts but some of them also disappear. 



3 For convenience we resort to the vocabulary pertaining to 3D triaxial tests 

4 When two grains arc in contact via 1, 2 or 3 contact points, only one contact is counted 



When a goes to 0.5, these transformations arc still active but with less intensity. Geometrical imbrications 
between clumps increase with a and are "more difficult to lose" during the biaxial tests. It is also interesting 
to observe that A seems to reach a threshold when a > 0.4, Fig. ??. This last observation can be correlated 
to the evolution of <f> p , which also reaches a threshold for the same value of a, Fig. ??. 

Focusing on A between the peak and the critical state, e a = e p and £& = e c , Fig. ??, we can observe 
that the increase of simple contacts is small for every a (A/g^n — 1 < 0.1) and complex contacts are mainly 
lost 1 — A/£i£ii > 0.1, especially when a is small. The greater a is, the smaller the proportion of complex 

contacts lost (grain imbrications are destroyed less) . This may be a clue that <pt of clumps increases with 
a, as seen in Fig. ??. Nevertheless, some new investigations on the evolution of contact orientations are 
proposed in the next section. 



6.2. Evolution of contact fabric for Cjj samples 

Contact direction and its evolution during the vertical compression tests are often analysed [351 ] . Focusing 
on the first part of the mechanical behaviour, from the isotropic state to the stress peak for example, it 
is well known that in dense samples, contacts are mainly lost in the extension direction and gained in the 
direction of compression, [36|. In figures ?? to ??, for two values of a, we present statistical analysis of the 
evolution of contact direction by the evaluation of V(9) = N Sb (9)/N ea (9), where e a and correspond to 
two successive vertical strains levels and where N Ex (9) is the number of contacts in the direction 9. V{9) = 1 
expresses that the number of contacts in the direction 9 remains constant between the two configurations 
studied. If V{9) < 1, contacts are lost and if V(9) > 1, contacts are gained in the direction 9. Integrated 
over 9, (V) is a global evaluation of the proportion of gained or lost contacts. 

We focus on the evolution of contact anisotropy between the isotropic state and the peak, Fig. ?? for 
a = 0.2 and Fig. ?? for a = 0.5, we can observe that there is no contact gain in any direction: in the 
compression direction, the number of contacts remains constant V(9) ~ 1 and the number of contacts 
decreases in the extension direction V(9) < 1. The mean value of V over 9 is smaller than 1 for both a 
discussed here and also for the other a studied. During the vertical compression of Crj , sample contacts 
are mainly lost. Finally, we noticed that the greater a is, the smaller the amount of contacts lost in the 
extension direction. By analysing the contacts change in direction between the peak and the critical state, 
e p to e c , Fig. ?? and ??, opposite tendencies emerge: contacts are mainly lost in the compression direction 
and gained in the extension direction, Fig. ??, with a less intensive effect for a = 0.5, Fig. ??. At this stage 
of analysis, we are not yet able to distinguish the nature of the contacts involved in these observations. An 
analysis of clustered contacts, as outlined below, is thus necessary. 

We now suggest the same contact direction analysis but for (SC) and (CC) groups (simple and complex 
contacts). Statistical analysis of the evolution of contact orientation from the isotropic state to the peak 
e p for (SC) and (CC) groups is shown in figures ?? to ??. On one hand, Fig. ?? and ?? show that in the 
compression direction (SC) contacts are gained (a = 0.2) or are kept (a = 0.5). (SC) contacts are mainly 
lost in the extension direction, with a more pronounced amplitude when a is small. On the other hand, we 
can observe that (CC), Figs. ?? and ??, are lost in every direction with some variations depending on 9. 
Nevertheless, complex contacts are more persistent when a is greater {{V) is greater for a = 0.5 because of 
grain imbrications). 

The statistical analysis of the evolution of contact direction between the peak e p and the critical state e c 
shown in the figures ?? to ?? confirms the tendency shown in figures ?? and ??: simple and complex contacts 
are lost in the compression direction. For (SC) with a = 0.2, (V) = 1, Fig. ??: although the number of 
(SC) decreases in the compression direction and increases in the extension direction, the number of simple 
contacts remains constant during the mechanical test from the peak to the critical state. When a = 0.5, the 
number of simple contacts decreases (V) = 0.9, especially in the compression direction. Complex contacts 
are the ones that are lost the most ((V) is always lower than 1, regardless of the value of a). Nevertheless, 
it is desirable to make a distinction based on a: the amount of (CC) lost is smaller for a = 0.5 ((V) = 0.8, 
Fig. ??) than for a = 0.2 ((P) = 0.7, Fig. ??). The proportion of (CC) lost in the compression direction is 
bigger for a = 0.2, like if vertical contact chains were more unstable or less persistent when a tends to 0. 



Further studies, taking into account the intensity of contact forces and their propensity to be stronger in 
the case of complex contacts (37j would provide more certainty about possible links between the observations 
made above and improvement of the mechanical property <p t measured at e c , Fig.??. 

6.3. Evolution of comparative contact proportions in compression tests for clumps and polygons 

Focusing on the critical phase only, contact observations are now based on the division into two contacts 
groups corner-to-edge (CE) and edge-to-edge (EE) (Fig. ??). The evolution of (CE) contacts according 
to a is shown in figure ??. (EE) contacts can be easily deduced by subtracting (CE) contacts percentage 
from 100%. On one hand, the percentage of (CE) contacts does not depend on the initial compacity of 
the sample; dense and loose samples exhibit a similar trend. On the other hand, ( CE) contact percentages 
for Cd and Cl samples decrease linearly with a. A different tendency is observed for Prj and Pl samples 
where the contact percentage remains more or less the same except for a — 0.5. We can even observe that 
when a = 0.5, (CE) contacts and (EE) contact percentages are close: for this special value of a, clump and 
polygon envelopes converge. 



6.4- Influence of shape on local strain analysis 

We focused on the strain localisation in the samples in order to study the macroscopic rupture and its 
origin. Two approaches were used: local strain maps and shear localisation indicator S2 [38[ . By comparing 
particle kinematics in the isotropic state E\ = £i = and in the deformed stage £1, we calculated local 
strains using Delauney triangulation as in (36^ (Delauney triangle corners correspond to the mass centers of 
particles). Using the second strain invariant, I £d = ej — en, where si and en are respectively the major 
and the minor principal strains, we illustrate the shear localisation in figure ??. We should notice that such 
shear localisation patterns, also called shear bands, are often observed on granular materials confined by 
more or less rigid boundaries, or even numerical or experimental considerations S [HEi Sl|. Even if 



the number of shear bands depends on the macroscopic strain levels applied to samples, different patterns 
(multiple shear zones) seem to exist when periodic boundary conditions are used [42j ■ 
The shear localisation indicator S% is defined as 




S2 = ^r v 'r. ' . (2) 



1=1 



E 7 " 



where l td is the second invariant of the strain tensor and Nt is the total number of Dclaunay triangles. In 
a sense, value of S2 can be regarded as a percentage of a distorted sample area. 

Figure ?? gives the evolution of S2 according to e\ for several dense and loose samples made of clumps 
and polygons. Regardless of the sample studied, the sheared area always reaches a maximum which is at least 
greater than 50%. From evolution of S 2 , we may encounter two types of behaviour: samples for which S2 
reaches a maximum and then decreases and stabilises, and a second type where S2 increases asymptotically 
towards a maximum. 

We have observed that when 5*2 reaches a maximum and later decreases to reach a threshold value, it 
always corresponds to samples which were identified as dense samples because <p p > <f>t (for example Cr>50 , 
Ppj.24 , Pl-28 samples of the figure ??). On the contrary, when 4> p ~ <f> t , samples can be classified as loose 
samples. In this case, S2 continuously increases from E\ = to 10% to eventually reach a threshold close to 
S 2 ~ 60%. 

In the case of dense samples, the asymptotic value of S2 is interesting because it clearly shows higher 
values for Prj samples than for Crj samples. Coupling this quantitative result with the qualitative observation 
of shear maps like in the figure ??, it is obvious that localised zones in Prj samples are wider than in 
Cd samples. This result is consistent with the overall dilatancy of samples: Prj samples globally expand 
more than Cd ones. For loose samples, the sheared area always corresponds to approximately 60% of the 
samples, regardless of the grain shape. 



7. Conclusions and discussion 

The aim of this article was to present some new investigations on the mechanical influences of particle 
shape in granular assemblies in the framework of numerical simulations performed with Discrete Element 
Method. First, a grain geometry parameter a was defined by the CEGEO research team. For particles 
called clumps, made of 3 overlapping discs, a is a measure of the grain concavity. 6-edged convex polygonal 
grains were also ruled by a. The overall envelope depending on a for each type of particles used in the 
studied granular model was the common feature. Our numerical simulations were performed with the 
Discrete Element Method adapted to each grain shape. For particles made of discs (clumps), the commercial 
code PFC 2D by ITASCA was used. For polygonal particles, we developed our own computer code which 
implements special contact detection between objects in the framework of Molecular Dynamic approach. 
In this article, we highlighted that changing the grain geometry, influences granular assembly mechanical 
behaviour under the classical vertical compression test in 2D. More complex grain shapes allow higher levels 
of internal friction angle and large volumetric strains to be reached compared to simple discs. Some clear 
differences in the behaviour of polygon (convex) and clump (non-convex) assemblies were shown. We should 
also note that the particle shapes chosen also demonstrate similarities, caused by the global envelope, that 
justify the comparison. The generation and compaction of granular assemblies was presented. By using 
two extreme values of the inter-granular friction angle fj,, dense and loose samples were prepared, both for 
samples made of clumps and polygons. 

Firstly, by focusing on the macroscopic mechanical behaviour of our granular model we show that loose 
samples composed of polygons with low values of a present behaviour typical to moderately dense granular 
samples. For these granular materials, the initial contracting stage was not only due to contact stiffness, 
but was also influenced by large inter-granular reorganisation. Apart from this, loose and dense samples 
of all shapes behaved as expected (loose samples only show contracting behaviour while dense ones mostly 
exhibit major dilatancy), showing similar behaviour when discussing friction residual angles <pt or contact 
percentages. All samples show higher values of internal friction angles <p p and </>t than samples made of only 
circular grains where each particle is a disc. The correlations between shape parameter a and friction angles 
are different for clumps and polygons. On one hand, for dense clump samples 4> p increases with a and seems 
to reach an asymptotic value <f> p = 40°. On the other hand, <p p linearly decreases when a shifts from 0.13 to 
0.5. In this case, the particular case of the triangular shape (a = 0.5) is also discussed briefly. The overall 
dilatancy of clump samples is greater than that of disc assemblies, but spectacularly smaller than dilatancy 
of polygons. 

Secondly, on the granular scale, we suggested correlating macroscopic observations by means of contact 
evolution analysis, which led us to introduce several groups of contacts between particles. Thus, we observed 
that multiple contacts between clumps transform into simple contacts and that this process depends on the 
size of concavities, i.e. a. We tend to associate this with an increase of shear resistance in the case of dense 
granular samples made of clumps. However, for a better understanding of the role of a in the magnitude 
of contact forces and their influence on the macroscopic repercussions, complementary studies need to be 
carried out. 

Focusing on granular assembly failure, we studied the localisation of shear bands and tried to characterise 
it by a scalar. It was observed that reflecting shear bands were thinner in dense samples made of clumps than 
in those made of polygons, regardless of the a type, thus highlighting the evident effects of the geometrical 
imbrications of clumps. Polygon samples gradually create wide shear bands, while for samples made of 
clumps, the appearance of shear bands are more immediate. 

Although the meaning and implications of the parameters a have been presented in this article, it needs 
to be clarified further. Nevertheless, from the comparisons between polygons and clumps, two trends seem to 
emerge: for very dense samples made of clumps, a large a naturally implies imbrications between particles. 
a is thus a measure of clump concavity. With the simulations exposed in this article, we can deduce that 
the larger the concavities are, the higher the angle of friction at the peak <j> p is. 

For the samples made of polygons, this does not apply. Indeed, there are no imbrications between the 
grains and a is thus a measure of grain sphericity. We have shown that in dense samples made of polygons, 
the percentage of a single contact in the isotropic state increases with a. In corollary, we can observe that 



the angles of friction peak decrease slightly with a. 

For loose samples PLOr Cl, or for Prjand Cd samples close to the critical state, fa is the parameter which 
characterises the failure. We have shown that for large strains, contacts between grains are mostly single. 
Therefore, clump imbrications are less involved in the evolution of fa. Furthermore, fa, which increases 
linearly with a, increases with the same rate regardless of the grain shape, a should thus be only regarded 
as a parameter of spherical grains. 
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